Rerouting Excitation Transfer in the Fenna-Matthews-Olson Complex 
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We investigate, using the Hierarchy method, the entanglement and the excitation transfer efh- 
ciency of the Fenna-Matthews-Olson complex under two different local modifications: the suppres- 
sion of transitions between particular sites and localized changes to the protein environment. We 
find that inhibiting the connection between the site-5 and site-6, or disconnecting site-5 from the 
complex completely, leads to an dramatic enhancement of the entanglement between site-6 and 
site-7. Similarly, the transfer efficiency actually increases if site-5 is disconnected from the com- 
plex entirely. We further show that if site-5 and site-7 are conjointly removed, the efficiency falls. 
This suggests that while not contributing to the transport efficiency in a normal complex, site-5 
introduces a redundant transport route in case of damage to site-7. Our results suggest an overall 
robustness of excitation energy transfer in the FMO complex under mutations, local defects, and 
other abnormal situations. 

PACS numbers: 



Photosynthesis is one of the most important bio- 
chemical processes on earth [l[. When light is absorbed 
by a light-harvesting antenna the excitation is trans- 
ferred to a reaction center and used for charge sepa- 
ration. Among the various photosynthetic complexes, 
the Fenna-Matthew-Olson (FMO) complex in green sul- 
fur bacteria is one of the most widely studied y. It 
has seven electronically-coupled chromophores and func- 
tionally connects a large light-harvesting antenna to the 
reaction center. Since the observation of quantum co- 
herent motion of an excitation within the FMO complex 
at 77 Kelvin [3|, considerable attention has been focused 
on the possible functional role of quantum coherence in 
photosynthesis [j|. Recent experiments further suggest 
the presence of quantum coherence even at room tem- 
perature |5[. 

Most quantum technologies, such as quantum compu- 
tation, quantum teleportation, quantum communication, 
rely on coherence in one way or another. Apart from 
photonic qubits almost all physical realizations demand 
extremely low-temperature environments to prevent fast 
dephasing [7[ and loss of quantum coherence. Therefore, 
the observation of quantum coherence (entanglement) in 
the FMO complex at ambient temperature has naturally 
triggered a great deal of theoretical interest and models 
[J, I8l-ll2| focusing on this biological system. The simplest 
theoretical treatment of the excitation transfer in the 
FMO complex normally considers seven mutually cou- 
pled sites (chromophores) and their interaction with the 



environment. One can either use the Lindblad master 
equation, the more accurate Hierarchy meth od J13l . or 
other open-quantum system models [9|, llO . Il4l |l5| to ex- 
plain the presence of quantum coherence and predict the 
physical quantities observed in experiments. 

In a natural in-vivo situation it is possible for the chro- 
mophores in the FMO complex to suffer damage, e.g., 
from optical bleaching or mutation, such that a transfer- 
ring pathway is blocked, or such that the environment 
(protein) is modified in some way. This has been demon- 
strated in recent experiments [l6j. Motivated by this 
fact, we investigate in this work how the entanglement 
and the transfer efficiency change when certain pathways 
are blocked, or the properties of the local environment of 
one site are modified. This question has been raised else- 



where, for example, Ref. [17[ discusses, using a Marko- 
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vian model, how various dissections of the FMO complex 
affect the efficiency and global entanglement. 

Here, we specifically focus on the situation where an 
excitation arrives at site-6, and must reach the reaction 
center at site-3 (similar roles may be played by site-1 and 
site-4, respectively). In this scenario, we ask the ques- 
tion what role is played by site-5 (see Fig. 1), and what 
happens if it, or site-7, are damaged? We find that if 
site-5 is damaged or removed from the complex entirely, 
the entanglement between sites 6 and 7 increases dra- 
matically, as docs the dynamic population of site-7 and 
consequently the efficiency (as characterized by the pop- 
ulation of the 'reaction centre') [18[. We then show that 
if site-7 is damaged conjointly with site-5, the efficiency 
falls. Thus, site-5, while not positively contributing to 
the efficiency in a perfect FMO complex, adds robustness 
and redundancy (as does the 6-1-2-3 transport route). 



We begin with an brief introduction to the standard 
model of the FMO complex, and the description of the 
environment using the Hierarchy equations of motion. 
We then discuss the concurrence and efficiency for dam- 
age and removal of site-5, and justify our interpretation 
of the role of site-5. Finally, we also consider a simplified 
Markovian model of a 3-site system and obtain analytical 
results for the concurrence between two of the sites, to 
further elucidate our full numerical data. 



Initial excitation 




FIG. 1: (Color online) Schematic diagram of a monomer of 
the FMO complex. The monomer consists of eight (only seven 
of them are presented here) chromophores. The excitation 
(from the light-harvesting antenna) arrives at sites 6 or 1 and 
then transfers from one chromophore to another. When the 
excitation arrives at site 3, it can irreversibly move to the 
reaction center. Here we assume that the initial excitation is 
at site 6. 



FMO MODEL 

Consider first a single FMO monomer containing N = 
7 sites, the general Hamiltonian of which can be written 
as 

N 

H = Y,^n\n){n\+ Y, Jn,n'i\n){n'\ + \n'){n\) (1) 



where the state \n) represents an excitation at site n 
(n G 1 ;...,?), Cn is the site energy of chromophore n, and 
Jn,n' is the excitonic coupling between the n-th and n'-th 



sites. Here, for simplicity, we omit the recently discovered 
eighth site [l9[ because its role on the excitation transfer 
process requires further studies. It has been shown that 
the exitonic coupling Jn,n' is of the same order as the 
reorganization energy, i.e., the coupling to the nuclear 
motion (phonons) of the protein environment. Thus a 
normal secular Redfield, or Markovian Lindblad, treat- 
ment is insufficient [111 l2l| , and the dynamics of the sys- 
tem must be modelled with a more complete approach, 
such as the Hierarchy equations of motion 13|. These 



equations are non-perturbative and non-Markovian, and 
valid under the assumption of a Drude spectral density 
and an initially separable system-bath state at t = 0. 
The Hierarchy is described by a set of coupled density 
matrices: 



N K \ N K 

Pn = ^ I -^ + X! X! nj.mMm Pn - i 

j—1 m—Q J j—1 m—0 

N K 
j—l m—0 






(2) 



Here, Qj = \j){j\ is the projector on the site j, L is 
the Liouvillian described by the Hamiltonian and the 
irreversible coupling to the reaction center (see below) 
L^ -jr[H, p„] + Lsink- Here, 



Cj = Y^ Cj,™ exp {-flj,mt) 



(3) 



rn=0 



where /ij.o = 7j, fJ-j.m > 1 = 27rr7i/ft/3, and the coeffi- 
cients 



',0 = Ij^j [cot{/3h'-fj/2) — i] /H 



and 
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(5) 



7j is the "Drude decay constant" , and indicates the mem- 
ory time of the bath for site j (each site is assumed to 
have its own independent bath), \j is the reorganisation 
energy, related to the system-bath coupling strength. 

A full description of the Hierarchy method can 
be found in the literature 13|, but in summary 
the Hierarchy is a large set of coupled equations 
each labelled by n, a set of non-negative inte- 
gers uniquely specifying each equation. The in- 
tegers are defined as n = {ni, n2, 71.3, ..., njv} = 
{{nio,nii,..,ni/f},..,{nAro,nAri,..,njvA:}}- In other 
words, each site j has an additional label m, from to 
K^ and each of those labels in turn can run from to 
00. The label n = = {{0, 0, 0....}} is special, and refers 
to the system density matrix. Its properties at any time 
t define those of the system. This is in turn coupled to 
so-called "auxiliary density matrices" , which describe the 



complex bath fluctuations, by the terms in the equation a simple thrcc-sitc model, which we discuss in section IV 
with n-^ (i.e., n -^ implies the term in the index dc- „ 
fined by j, and m is increased or decreased by 1). At 
high temperature, and imposing the Ishizaki-Tanimura 
boundary condition [13| , we can cut the Hierarchy off at 
K = and an appropriate total number of terms in the 
remaining labels Nc = J2j m '"'i,™ providing convergence. 
We also include isink to describe the irreversible exci- 
tation transfer from site 3 to the reaction center: 



Lsink[p] = r[spst 






(6) 



Efficiency in the presence of defects 

What do these concurrence results imply for the overall 
efficiency of the transport process? The efficiency of the 
excitation transfer can be formulated via the population 
of the reaction center as a function of time, 



where s = |0)(3|, with |0) denoting the state of the reac- 
tion center, and F the transfer rate. 

In the FMO monomer, the excitation transfer from site 
3 to the reaction center occurs on a time scale of ~ 1 ps, 
and the dephasing occurs on a time scale of ^ 100 fs [12| . 
This two time scales are both much faster than that of 
the excitonic fluorescence relaxation (^ 1 ns), which, for 
simplicity, is omitted in our explicit results. 



CONCURRENCE AND POPULATION 
DYNAMICS IN THE PRESENCE OF DEFECTS 

Each site in the FMO monomer may be decoupled from 
its nearest-neighbor sites due to mutation-induced de- 
fects or rotation of the site [l6[ . To investigate the effect 
of this change on the excitation transfer we consider the 
situation where the initial excitation arrives at site-6 and 
study the temporal excitation transfer when the excitonic 
coupling between two specific sites is inhibited. Wc find 
that when the coupling between site-5 and site-6 is inhib- 
ited, a significant enhancement of the coherence between 
site-6 and site-7 can be obtained. To examine the co- 
herence between the two sites, we utilize the bipartite 
concurrence C, which quantifies the degree of entangle- 
ment of any two sites n, n': 



Cn.n'=2\{n\po\n')\. 



(7) 



This is extracted from the po density matrix, evaluated 
from the time evolution of the Hierarchical equations of 
motion. 

In Fig. 2(a,b), we show the concurrence C of site-6 and 
site-7 when only the excitonic coupling between site-1 
and site-6 is inhibited [blue dashed line, for bath temper- 
atures of 77 K in (a) and 300 K in (b)]. The concurrence 
increases slightly since less sites share the excitation from 
site-6. In contrast, when the coupling between the just 
site-5 and site-6 is inhibited, a much larger enhancement 
(red dashed curve) of the coherence can be obtained. 
This is simply because the 5-6 coupling is much larger 
than the 6-1 coupling, thus when it is inhibited more 
population can flow to site 7, which, since it is a co- 
herent process even at 300 K, increases the concurrence 
between sites 6 and 7. This can be further clarified with 



Pnc{t)=Tr[p{t)ss^ 



(8) 



where s = |0)(3| is the operator connecting site-3 to the 
reaction center (see Fig. 1), denoted by the state |0), as 
defined earlier for the Lindblad Lsink- Since the excitonic 
fluorescence relaxation of each individual site is slow {^ 1 
ns) compared to all other time scales, Pncit -^ oo) ap- 
proaches unity, leading to the near 99% efficiency of the 
FMO complex commonly discussed in the literature. It 
is often argued that coherence plays an important role 
promoting this high-efficiency, but some interesting in- 
vestigations have shown that the quantum and classical 
models only differ by a few percent [23] ■ 

To check this long-time behavior we employed an ex- 
tended model (results not shown here), including the ex- 
citonic recombination rate of each individual site, and 
found that the defects discussed in the previous section 
do have a small effect on the long-time dynamics, and 
that the magnitude of this effect strongly depends on 
parameters which are not precisely known [8|. For ex- 
ample, a change in the rate F between site-3 and the 
reaction center by a factor of 5 results in a magnifica- 
tion of any differences in the efficiency. In addition, we 
found that any such change in the long-time dynamics is 
pre-determined by larger changes in the early-time pop- 
ulation of the reaction center. Thus, here we use these 
short-time dynamics, in the absence of excitonic recom- 
bination, as an indicator of the efficiency. 

In figure 3 we show the reaction center population as 
a function of time for a range of defects. We see that, 
at both 77 K and 300 K, completely cutting site-5 (red 
dashed curve) enhances the reaction center population 
over the unmodified case (solid black curve), and hence 
enhances the efficiency. Conversely, cutting site-7 alone 
(long dashed blue curve) reduces the efficiency; in this 
case the population is forced to traverse through site-5, 
which is a less efficient, and slower, route to site 3, and 
subsequently to the reaction center. In contrast, remov- 
ing both site-5 and site-7 (green dotted curve) leaves only 
the 6-1-2-3 transfer route, which due to the weak coupling 
between site 6 and 1 is less efficient. This supports our 
earlier hypothesis that while site 5 does not contribute 
in a positive way to a perfect FMO complex, it does pro- 
vide necessary redundancy in case of damage to the more 
efficient transport through site 7. 
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FIG. 2: (Color online) (a,b) The concurrence (coherence) of sites 6 and 7 for (a) a bath temperature of 77 K and (b) 300 K. 
Not surprisingly, higher temperatures suppress oscillations. When the excitonic coupling between site-1 and site-6 (blue dashed 
curve) is inhibited the concurrence increases slightly. When the coupling between between site-5 and site-6 is inhibited (red 
dashed curve) or when site-5 is completely removed from the complex (orange dotted curve) the concurrence between 6 and 
7 rises drastically. The solid black curve represents the concurrence of the site-6 and site-7 for the full unmodified complex. 
Figures (c,d) show the behavior of the populations at 300 K. Figure (c) shows the full unmodified dynamics, while (d) shows 
the case where site-5 is completely removed, and hence the population of 7 rises at a faster rate. Interestingly, in (d) the 
coherent oscillations in the site-6 population disappear, while in (b) we see that the concurrence remains large, indicating that 
sometimes coherent oscillations are not a strong indicator of coherence (as also seen in Q^). In plotting this figure, we set 



50 fs and A 



35 cm , and the rate from the site 3 to the reaction center F 



1 ps. 



LOCAL ENVIRONMENTAL CHANGES 

Other localized phenomena can also affect the transfer 
kinetics of photosynthetic complexes. For example, the 
vibronic structure of the FMO complex can be altered 
by the local substitution or deletion of the gene-encoding 
enzyme responsible for reducing the isoprenoid tail of the 
chromophores jl6|. These alterations can lead to modifi- 
cations both of the protein and the chromophores. 

We can easily investigate the effect on the coherence 
when the local environment of one site is changed. We 



assume that the modification of the local environment of 
site-5 results in a stronger coupling to the protein envi- 
ronment. In Fig. 4, we show the concurrence C of site- 
6 and site-7 as the coupling site-5 to the environment, 
via the reorganisation energy A, is increased. In con- 
trast to when site-5 was removed from the complex, the 
concurrence initially increases slightly, due to a small en- 
hancement of the population flowing to site-7, and then 
decreases (with respect to the unmodified complex). We 
also observed the change in the reaction center popula- 
tion as a function of these changes in A (not shown in 
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FIG. 3: (Color online) (a,b) The population of the reaction center (an indicator of the transport efficiency) for (a) a bath 
temperature of 77 K and (b) 300 K. When site-5 is removed from the complex (red dashed curve) the efficiency actually rises. 
When site-7 is removed (blue dashed curve), the population falls, while when both 5 and 7 are removed the population falls 
even further (green dotted curve). The solid black curve represents the population of the reaction center for the full unmodified 



complex. As before, we set 7 



50 fs and A 



35 cm , and the rate from the site 3 to the reaction center F 



1 ps. 



the figure) and found an overall small increase, but not 
as drastic as that observed in Fig. 2. The dynamics and 
efficiency, as a function of global changes in the environ- 
mental coupling and damping rate, have been well stud- 
ied in many other works |lll . |22| . A full investigation, 
more rigorously taking into account the physical effect of 
the changes observed in [l6|, remains to be performed. 
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FIG. 4: (Color online) Concurrence of site-6 and site-7 for 
increased coupling of site-5 to the protein environment, A5 — 
35 cm^^ (black solid), A5 = 70 cm~^ (red dashed), and A — 
175 cm~ (green dashed). The bath temperature of all sites 
is kept at 300 K, and 7"^ = 50 fs, F~^ = 1 ps, as before. 



THREE-SITE MODEL 

To further understand the mechanism that leads to the 
enhancement of the coherence between sitc-6 and site-7 
observed in Fig. 2, we now consider a simplified three- 
site model, as shown in Fig. 5, where the three sites, 1, 
2, and 3, now represent site-6, -7, and -5, respectively, 
in the FMO complex. We assume that the initial ex- 
citation is at site 1, and assign the inter-site couplings 
values approximately corresponding to the excitonic cou- 
plings in the FMO monomer: Ji = Je^r, J2 = ^5.7, and 
J3 = J^fi. We further apply a Markovian dissipative 
channel to both sites 2 and 3 with the rate 7 to simu- 
late the excitation transferring to other sites in the FMO 
monomer. Because J2 is small compared with Ji and J3, 
we can approximately set J2 = 0. For simplicity, we fur- 
ther assume Ji = J and J3 = ^J, where ^ is the tuning 
parameter for the coupling strength. The concurrence of 
this simplified model can then be expressed as 



C 



where 



je-*(7+r)(ert_i)j 



2fir3 



+ 7[e^*(7 + r)]} 



{F^ + (f2 - 7^)6 



rt 



v/-4(l + e2)(J/ft)2+72. 



(9) 



(10) 



As shown in Fig. 5(a), the concurrence of sites 1 and 
2 is strongly enhanced (black solid curve) when the cou- 
pling J3 between sites 1 and 3 is switched off. This is be- 
cause the excitation population is predominantly trapped 



between sites 1 and 2. In Fig. 5(b), we further show how 
the concurrence increases while decreasing the couphng 




FIG. 5: (Color online) A simple model of three sites simulates 
the enhancement of the coherence observed in the full model 
of the FMO complex in Fig. 2(a, b). Here the sites 1,2, and 
3 represent sites-6, -7, and -5 in the FMO monomer, respec- 
tively. The model Hamiltonian has the same form (though 
now three sites are included) as Eq. (1), and the inter-site 
couplings are set to be Ji = Jg,?, J2 = Js.r, and J3 = Js.e. 
Because J2 is small compared with Ji and J3, we approxi- 
mately set J2 = 0. For simplicity, we further assume ,J\ = J 
and Ji = ^J with ^ being the tuning parameter. In order 
to simulate the excitation transfer simply, we apply an addi- 
tional Markovian dissipative channel with rate 7 to both sites 
2 and 3. In (a) the red curve shows the concurrence between 
sites 1 and 2 for the full 3-site model. The black curve shows 
the concurrence when site-3 is decoupled from site-1. Figure 
(b) shows the variation of the concurrence between site-1 and 
2 as a function of time and the parameter ^ of J3 . In plotting 
this figure, we apply 7 — 5.3 cm^^ and J — 30cm^^ . 



CONCLUSIONS 

In summary, we have shown that the concurrence (co- 
herence) between site-6 and -7 can be enhanced signif- 
icantly when the coupling between the sites-5 and -6 is 
inhibited, or if site-5 is removed from the complex com- 
pletely. In the latter case we also found a corresponding 
increase in the population of the reaction center. We 
then argued that, rather than being superflous (and in 
fact contributing negatively to the efficiency of a perfect 



FMO monomer), site-5 provides a backup in case of dam- 
age to the highly-efficient transport through site-7. We 
further apply a simplified three-site model to simulate 
this result. Overall our results imply a robustness and 
redundancy to the energy transfer in the FMO complex, 
as also noted in jl7|. 
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